Theory for the ultrafast ablation of graphite films 
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The physical mechanisms for damage formation in graphite films induced by femtosecond laser 
pulses are analyzed using a microscopic electronic theory. We describe the nonequilibrium dy- 
namics of electrons and lattice by performing molecular dynamics simulations on time-dependent 
potential energy surfaces. We show that graphite has the unique property of exhibiting two dis- 
tinct laser induced structural instabilities. For high absorbed energies (> 3.3 eV/atom) we find 
nonequilibrium melting followed by fast evaporation. For low intensities above the damage thresh- 
old (> 2.0 eV/atom) ablation occurs via removal of intact graphite sheets. 



in 



52.38. Mf, 64.70.Hz, 



;.35.Ja, 81.05.Uw 



C/2 



o 



> 
00 



o 



I 

O 

o 



Material irradiation with femtosecond laser pulses 
gives rise to many interesting phenomena involving 
atomic motion, like ultrafast phase transitions ab- 
lation |^,|| , optical breakdown |p^|~|l^ , and excitation of 
coherent phonons . Ultrafast laser ablation, in partic- 
ular, has been extensively studied mainly because of its 
important technological applications [Q. However, little 
is known about the microscopic mechanisms leading to 
ablation. 

We define the ablation threshold as the laser fluence for 
which lattice instabilities of such magnitude are induced 
that the system is irreversibly damaged and at least a 
monolayer of material is removed [|o|. On the basis of 
many experimental studies it is commonly assumed that 
the ablation process near the threshold is always initiated 
by the melting (or ultrafast melting) of the material. 

In this letter we show that there is an exception to this 
general rule. We present a theoretical study of the mi- 
croscopic processes leading to the femtosecond ablation 
of graphite films. We conclude that, due to its layered, 
sp^-bonded structure, graphite presents the unique prop- 
erty of exhibiting two different ablation mechanisms and 
therefore two different ablation thresholds. The mecha- 
nism for the low fluence ablation threshold is the removal 
of intact graphite sheets and does not involve melting. In 
contrast, the high fluence threshold corresponds to bond 
breaking processes inside the graphite layers and leads 
to ultrafast melting and expansion of the film. More- 
over, we find that the metastable and short-lived low den- 
sity liquid state formed at the second ablation threshold 
presents some signatures corresponding to low density 
liquid carbon (LDLC) |l|l. 

We study the laser-induced ultrafast dynamics of a 
graphite film by performing molecular dynamics calcu- 
lations based upon a tight-binding Hamiltonian. The ef- 
fects of the strong electronic nonequilibrium created by 
ultrashort laser pulses are accounted for with a method of 
determining time-dependent electronic occupation num- 
bers, which leads to potential energy surfaces changing 
with time p]. The film geometry is simulated with a 



MD supercell that has periodic boundary conditions in 
the horizontal directions, while in the direction perpen- 
dicular to this, the material borders on the vacuum [ p^ . 
The shape of the MD supercell is constant in this model, 
and the size is sufficiently large to allow any expansion 
of the material along the vertical axis. With this method 
we treat the inertial confinement which is expected to 
be important in the initial stages of the response of a 
graphite film to laser excitation. The method of using 
a column-shaped MD supercell of fixed shape and size 
should be a good approximation to the physical reality 
in the first 10 to 100 ps after laser excitation, i. e. as long 
as the vertical expansion of the material dominates and 
the horizontal expansion can still be neglected. 

In the MD simulations the forces acting on atom k is 
calculated as ik{{rij{t)},t) = — Vfe $({ry (i)}, t), where 
rij = \ri — Tjl is the distance between atoms i and j, 
and ^{{rij{t)},t) is the potential. It is clear that, for 
the description of laser induced nonequilibrium processes, 
$({ry (t)},t) cannot be a model potential, but it has to 
be derived from a microscopic electronic theory, and is 
given by 

'^i{nj{t)},t) = - ^n(e„,t)e™ - ^ ^ £;rcp(r.j) . (1) 

m ij 

In Eq. (|l|) the quantities e™ — {m\H'TB{{'''ij{t)})\m) are 
the eigenvalues of the electronic Hamiltonian, n{em,t) 
the corresponding electronic occupations, and i?rop(''ij) 
the repulsive core-core potential. We use a tight-binding 
Hamiltonian 
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Here, represents the occupation number operator for 
the orbital 77 of atom i, and the hopping matrix ele- 
ments have been abridged by tj^ . We follow Slater and 
Koster [|l^ in their treatment of the direction dependence 



We consider the four valence 



of the hopping integrals t/^ 
orbitals of carbon 2s, 2px, '^Py and 2pz (thus rj^C. — ssa, 
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spa, ppa and ppn). For the distance dependence of the 
TB parameters Vssa, Vspa, Vppcr and YppT^, and the re- 
pulsive potential the form proposed by Ho et al. |]l8| is 
adopted. 

Our treatment of the effect of the optically created 
electron- hole plasma is founded on the following physical 
picture: Due to the action of the laser pulse electrons are 
excited from occupied to unoccupied levels with a time- 
dependent probability which is proportional to the inten- 
sity of the laser field. As a consequence of the extremely 
fast excitation process, a non-equilibrium distribution of 
electrons is created. Through electron-electron collisions 
this electron distribution thermalizes to an equilibrium 
occupation of the electronic levels. We describe the ab- 
sorption of energy by the electronic system and its equi- 
libration by: 

+n{em + hcj,t- At) - 2n{em, t - At)] | 

n 

(3) 

Thus, the electronic distribution is at each time step 
folded with the pulse intensity function g{u},t). This 
means that at each time step, the occupation of an energy 
level em changes in proportion to the occupation differ- 
ence with respect to levels at e,„ — Klo and at e„i + huj. 
In Eq. (^, constant optical matrix elements are assumed. 
The second term of Eq. (||) describes the electron-electron 
collisions that lead to an equilibration of the electronic 
system with a rate equation of the Boltzmann type for 
the distribution n(em, t). Hence, with a time constant ti, 
the distribution n{em, t) approaches a Fermi-Dirac distri- 
bution nP{€„L,Tc) (with electron temperature To). 

This simple approach works very well because for dense 
electron- hole plasmas extremely short relaxation times ti 
have been found. Chemla and coworkers have reported 
a carrier thermalization faster than 10 fs in GaAs [p^ . 
As we are not aware of a measured relaxation time in 
graphite, we use ri = 10 fs. Note, for such a short ther- 
malization time the exact electronic dynamics leading 
to electronic equilibrium do not play a significant role 
for the structural changes we are studying here. The 
electronic temperature To and the chemical potential /i, 
which appear in the Fermi-Dirac distribution and which 
are not determined by Eq. (|^), need to be fixed by an ad- 
ditional principle. We demand that the noncquilibrium 
distribution n(e,„,i) approaches the Fermi-Dirac distri- 
bution while conserving the total energy of the system. 
On the extremely short time scale of ti = 10 fs, we expect 
energy loss mechanisms to be negligible. 



Using the theory described above we perform simula- 
tions on graphite films for a wide range of pulse intensi- 
ties and durations. This allows us to determine ablation 
thresholds. 

In Fig. we show snapshots of a long trajectory of 
a graphite film in which melting and subsequent evapo- 
ration can be observed. The film consists of = 576 
atoms. The dynamics shown in this and the following 
figures are for laser pulses of t = 20 fs duration with a 
central frequency of Tiuj = 3.0 eV. We also calculated tra- 
jectories at fiuj = 2.0 eV and found the same dynamics. 
The intensity was determined so that the film absorbs 
Eq = 4.0 eV/atom, corresponding to a carrier density of 
= 6.1 X 10^^ cm~'^ and a fiuence of F « 0.35 J/cm^ 
|pO| . This energy density is high enough to break the 
bonds of the graphite planes. At t = 80 fs after the 
pulse maximum (see Fig. ^ (b)) we see that disorder has 
developed inside the graphite planes and first bonds are 
formed between them. At the same time, the first car- 
bon monomer is evaporated from the material, confirm- 
ing that the graphite planes have been excited above their 
fragmentation threshold. The following snapshots show 
the complete disintegration of the graphite planes. In 
Fig. (c) to (e), a strong volume expansion is observed. 
In Fig. |l| (c) to (e), four monomers that had already 
moved far away from the surface were omitted. In the 
advancing ablation process, further carbon monomers are 
ejected, and between At — 520 fs and At = 920 fs we can 
observe how the expansion of the material takes place 
especially rapidly in the two surface regions of the film, 
leading to low densities and the formation of short car- 
bon chains which then start to leave the sample. For the 
ablation of graphite via the mechanism described above 
we determined a threshold of tab2 = 3.3 ± 0.3 eV/atom, 
corresponding a fiuence of i^Q;,2 ~ 0.29 ±0.03 J/cm^ [^ . 

The examination of a graphite film which had absorbed 
a moderate energy density has produced a surprising re- 
sult: There is an ablation mechanism which does not lead 
to a destruction of the graphite planes. The physics of 
this mechanism is as follows: First, in thermal equilib- 
rium the graphite planes are separated by a very large 
distance of c? = 3.4 A and consequently their interaction 
consists only of a very weak van-der-Waals-interaction 
energy E Ki \2 meV Then a photoinduced electron- 
hole plasma causes a very strong vibrational excitation of 
the graphite layers and at the turning point of their oscil- 
lations the atoms of two different layers have a distance 
of d « 2 A. At this distance the interaction becomes al- 
ready quite strong and is on the order oi E m 1 eV. Thus, 
the laser-induced strong vibrations can lead to collisions 
of the planes in which momentum in the z direction is 
transferred. As a result, a surface plane of a graphite 
sample which usually has zero total momentum can gain 
enough momentum to leave the surface. 
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FIG. 1. Ablation of a 
graphite film for an absorbed 
energy of Eo = 4.0 eV/atom. 
The laser pulse duration was 
T — 20 fs. Note the strong 
expansion, the formation of 
a metastable liquid-like state 
and the emission of carbon 
atoms and chains. 



FIG. 2. Ablation of 

graphite for an absorbed en- 
ergy of Eo = 2.4 eV/atom, 
below the threshold for the 
destruction of graphite 
planes. The pulse duration 
was r = 20 fs. Thus, the laser 
pulse induces strong vibra- 
tions of the graphite planes 
which lead to collisions of the 
planes. As a consequence, 
the planes at the top and at 
the bottom are removed from 
the surface of the film. 



In Fig. H we present snapshots of the dynamics of 
graphite for a r = 20 fs laser pulse with an absorbed 
energy of Eq = 2.4 eV/atom, corresponding to a car- 
rier density of ric — 3.5 x 10^^ cni~^ and a fluence of 
F « 0.21 J/cm^ Q. The number of atoms is = 576. 
This trajectory shows the new ablation mechanism. At a 
deposited energy density that is below the fragmentation 



threshold of the graphite planes tab2, we find the abla- 
tion of entire graphite planes by the mechanism explained 
above: The laser pulse leads to a strong vibrational exci- 
tation of the graphite layers, with a significant movement 
of the atoms perpendicular to the graphite planes. This 
excitation is already strong at f = 80 fs after the laser 
pulse maximum, but for the topmost plane it reaches a 
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maximum at t = 160 fs (see Fig. ||). At that time the 
distance between the graphite planes has decreased to 
approximately d = 1.8 A, and the first and second planes 
strongly interact. We can speak of a collision of the two 
topmost graphite planes, which then leads to a momen- 
tum transfer. Consequently, the surface plane that had 
up to that time zero total momentum starts to leave the 
crystal in positive z direction. In the same way the bot- 
tom layer of Fig. ^ can be seen to leave the rest of the 
film in negative z direction. 

This ablation mechanism in graphite has not been de- 
scribed before. It solely relies on the fact that strong 
vibrational excitation of the graphite planes can increase 
the interaction between the planes to such a degree that 
the momentum transferred in this collision process causes 
ablation of entire planes. It is to be expected, of course, 
that in reality we are dealing with graphite planes that 
show slight defects. These would result in large plane 
fragments being removed from the surface instead of 
entire graphite planes as the simulation with periodic 
boundary conditions in the horizontal directions sug- 
gests. We estimate the threshold for this ablation mech- 
anism to be tab = 2.0 ± 0.4 eV/atom |E2|, corresponding 
a fluence of Fat « 0.17 ± 0.04 J/cm^ [g. 

The absorbed energy Eq — 2.4 eV/atom in the exam- 
ple shown in Fig. || is high enough to cause the ejection of 
the first graphite sheet within one picosecond. Note, this 
absorbed energy is large enough to produce a thermal 
melting of the removed graphite planes at longer times. 
This means that although the ablation does not occur 
via melting, the ablation products might undergo melt- 
ing after some time. 
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FIG. 3. Pair correlation function g{r) during the 
melting of a graphite film for an absorbed energy of 
Eq — 4.0 eV/atom. The pulse duration was r = 20 fs. The 
pair correlation functions correspond to the same trajectory 
as Fig. |l| Note, the ordinate scale of (a) is different from that 
of (b) to (d). 



We now analyze in detail the metastable liquid phase 
induced above tab2, corresponding to the snapshots of 
Fig. |l|. In Fig. || we show the pair correlation functions 
g{r) during the ultrafast melting of a graphite film. They 
are calculated as 



lir) 



Ua 47r Ar 



(4) 



where AN{r) is the number of atoms situated at a dis- 
tance between r and r + Ar around a given particle, and 
Ua = N/il is the atomic density of the material. g{r) 
at the pulse maximum (see Fig. || (a)) docs not yet dif- 
fer much from g{r) before the pulse. But then a rapid 
broadening of the structures of g{r) occurs (Fig. ^ (b)) 
and at < = 80 fs the pair correlation function of liquid 
carbon is recognizable (Fig. ^ (c)). Note, because of the 
expansion of the film this liquid phase has a low density 
and is characterized by the presence of carbon chains. 
To characterize this laser induced liquid phase we have 
determined the order parameter 5* = (p4 — p2)/{pA + P2) 



which was proposed by Ree et al. |15| to distinguish be- 
tween the high density and the low density liquid phase 
of carbon. Here, p2 and p4 correspond to the densities of 
2-fold and 4-fold coordinated carbon atoms determined 
following Brenner |^ . At a time t = \ ps after the pulse 
maximum, we find a value of 5* = —0.8 for the graphite 
film shown in Fig. Q confirming that a low density, pre- 
dominantly 2-fold coordinated liquid carbon phase was 
produced by the laser excitation. 

For a clear check of our predictions pump-probe ex- 
periments with a time resolution of less than 100 fs and 
a structural resolution of a few atomic layers would be 
required. Such experiments have not been performed 
on graphite so far. However, the time-resolved reflectiv- 
ity measurements of laser-irradiated graphite that were 
reported by Sokolowski-Tinten et al. could be un- 
derstood qualitatively with the help of our calculations. 
As in the experiment we can distinguish two different 
ablation thresholds {tab,tab2)- Our calculated ablation 
thresholds arc in good agreement with the experimental 
ones Fab = 0.185 J/cm^ and Fab2 = 0.25 J/cm^ [||. 
This demonstrates the validity of our approach using a 
tight-binding Hamiltonian. It is important to point out 
that none of the existing ah initio approaches to laser 
induced melting are able to take the pulse duration and 
the nonequilibrium electron dynamics into account. Note 
also that the number of = 576 atoms in the MD su- 
percell is not accessible to ah initio calculations PJl^]. 

An important result of our approach is that both ab- 
lation thresholds are nearly independent of the pulse du- 
ration for a range [10fs-500fs]. For fixed pulse duration 
and increasing energy, qualitative changes only occur at 
the ablation thresholds. Both the emission of graphite 
sheets and the ultrafast melting are initiated by strong 
in- and out-of-plane oscillations of the atoms. 
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Summarizing, the microscopic processes leading to 
laser ablation in graphite films have been studied the- 
oretically. A new ablation mechanism has been observed 
at energies below the threshold for the destruction of 
graphite planes. This ablation process is very different 
from ablation at higher deposited energies, since large 
graphite plane segments instead of small carbon clus- 
ters (chains) are produced. The strong anisotropy of 
the graphite structure is crucial for this ablation without 
melting. This process may play a role in the formation 
of carbon nanotubes by graphite ablation. 
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